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ABSTRACT 

Using K and R band imaging of the ELAIS N2 field, we investigate the number counts, 
clustering, morphology and radio/X-ray emission of extremely red objects (EROs), 
[ defined as galaxies with R — K > 5.0. This criterion will select old, passive ellipticals 

■ at z > 0.9. To K = 21 we identify a total of 158 EROs in 81.5 arcmin 2 . The ERO 
Q\ ' number counts are lower than predicted by pure luminosity evolution models, but 
I/"") , higher than predicted by current CDM-based hierarchical models. The ERO counts 
CNl ■ are consistent with a non-evolving model and also with a luminosity evolution model 

incorporating moderate merging and a decrease with redshift in the comoving number 
density of passive galaxies ('M-DE'). 

We investigate the clustering of the EROs by calculating their angular correlation 
function, u>(8), and obtain a > 2a detection of clustering at K = 19-20 limits. The 

■ a; (6*) amplitude of these EROs is much higher than that of full iif-limited samples 
Qu of galaxies, and best-fitted by models with a comoving correlation radius ro ~ 10- 

13 h~ x Mpc. These results, which are in agreement with Daddi et al. (2000), suggest 
£^ ■ that the intrinsic clustering of at least the brighter EROs is even stronger than that 

^ \ of present-day giant ellipticals. 

C$ We estimate seeing-corrected angular sizes and morphological types for a K < 19.5 

subsample of EROs (31 galaxies) and find a ^3:2 mixture of bulge and disk profiles. 
Of these EROs ~ j appear to be interacting, disturbed or otherwise irregular, and two 
>^ . are visible mergers. We find the angular sizes of the bulge-profile EROs are consistent 

$_i ' with passively evolving ellipticals in the M-DE model, at the expected z ~ 0.9-2. 

The ERO mean radii are smaller than the non-evolving prediction, implying surface 
brightness evolution. 

Seven of the 31 bright EROs are detected as F{XA GHz) > 30/iJy radio sources 
in a VLA survey. The strongest, at 5 mJy, is also a Chandra X-ray detection, and lies 
at the centre of a significant overdensity of EROs - it is probably an FRI radio galaxy 
in a z ~ 1 cluster. Of the other, much fainter, sources, five are point-like and may 
be weak AGN, while the sixth is elongated and aligned with the optical axis of an 
extended, low-surface brightness ERO, and hence probably a ULIRG-type starburst. 

A possible interpretation is discussed in which the EROs are a mixture of (i) 
'pEROs', strongly clustered passively evolving giant ellipticals, formed at high red- 
shifts, the oldest EROs, and (ii) 'dsfEROs', dusty post-interaction galaxies, with a 
few active starbursts (ULIRGs), and less strongly clustered. With time, the younger 
dsfEROs are continually assimilated into to the ERO class, diluting the clustering and 
increasing the comoving number density. Both types ultimately evolve into today's 
early- type galaxies. 

Key words: galaxies: evolution; galaxies: elliptical and lenticular, cD; galaxies: high- 
redshift; radio continuum: galaxies 
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1 INTRODUCTION 

Deep surveys at near-infra-red (~ 2/j.m) wavelengths re- 
vealed that a population of very red (R—K > 5 or I—K > 4, 
Vega system) galaxies appears faintward of K ~ 18 (e.g. El- 
ston et al. 1988; Hu and Ridgeway 1994). The colours and 
K magnitude range of these 'extremely red objects' (EROs) 
are consistent with old, passive (i.e. no longer forming new 
stars) galaxies, observed at z ~ 1 or beyond. These would 
be the progenitors of present-day giant ellipticals, thought 
to have formed their stars at even higher (z > 3) redshifts. 

The EROs are of great cosmological interest, as they 
will include the earliest formed of all galaxies, and their 
properties will help to answer the long-standing myster- 
ies concerning the early evolution of elliptical galaxies (e.g. 
Jimenez et al. 1999). One hypothesis is that the present-day 
ellipticals formed in single starbursts at z > 3, thereafter un- 
dergoing only passive luminosity evolution, as in pure lumi- 
nosity evolution (PLE) models, another is that they formed 
more recently through mergers of spiral galaxies, or by some 
combination of the two. Initial surveys of EROs (e.g. Barger 
et al. 1999) found fewer than expected from PLE models, 
and implied that < 50 per cent of the local E/SOs could 
have formed in a single high-redshift starburst. However, 
the detection with SCUBA of many sub-mm sources, with 
optical magnitudes and colours consistent with dusty galax- 
ies at z > 2, and 850^im fluxes indicative of star-formation 
rates (SFRs) as high as ~ lOOOM yr _1 , provided strong ev- 
idence that some massive galaxies did form in intense early 
starbursts (e.g. Smail et al. 1999). The passive EROs could 
then be the intermediate stage between the SCUBA sources 
and local giant ellipticals. 

With sufficient dust-reddening, high-redshift star- 
forming galaxies may also have R — K > 5. Spectroscopy 
of a few of the brightest EROs revealed that both old, pas- 
sive (e.g. Dunlop et al. 1996; Spinrad et al. 1997; Stanford et 
al. 1997) and dusty starburst (e.g. Dey et al. 1999; Smith et 
al. 2001) galaxies were present, but the relative proportions 
remained uncertain. 

The morphologies of the EROs provide further informa- 
tion on their nature. Moriondo, Cimatti and Daddi (2000) 
fitted radial profiles to WFPC2 and NICMOS images of 41 
EROs (to K ~ 21) and estimated that 50-80 percent were 
elliptical-like and 15 percent irregular or interacting. Sti- 
avelli and Treu (2000) studied 30 EROs to H ~ 23, using 
NICMOS, and classified these into E/S0, disk, irregular and 
point-like morphologies in a ratio 18:6:3:3. They assumed 
that only the first type were the progenitors of local E/SOs, 
and by comparing the counts of these with a PLE model 
(with Q m — 0.3, Qa = 0.7), estimated only ~ 15 per cent 
of the present-day number of E/SOs could have formed at 
z > 3. 

Daddi et al. (2000) presented the first measurements of 
ERO clustering. They found the angular correlation func- 
tion, uj{9), of K < 19.2 EROs to be almost an order of 
magnitude higher than that of all galaxies to the same K 
limit. As local giant ellipticals are intrinsically much more 
clustered than disk galaxies (e.g. Guzzo et al. 1997), this was 
interpreted as evidence that most EROs are z > 1 ellipticals. 

Manucci et al. (2001) attempted a photometric sepa- 
ration of passive and dusty, star-forming EROs - hereafter 
pEROs and dsfEROs - on a plot of J — K against R — K. 



Within a sample of 57 EROs with K < 20, they assigned 
equal numbers (21) to each class. Cimatti et al. (2002a) pro- 
vided the first true measurement of the dsfERO/pERO ra- 
tio, identifying by spectroscopy approximately equal num- 
bers (15 'dsf, 14 'p') amongst 45 EROs with K < 19.2. 
Firth et al. (2002) investigated ERO clustering to H = 20.5 
(K ~ 19.5) and obtained similar results to Daddi et al. 
(2000). 

In this paper we select a sample of R — K > 5.0 EROs 
from new K and R band images, investigate their number 
counts, clustering and morphologies and compare these re- 
sults with evolutionary models. We do not yet have spec- 
troscopy for these galaxies, but do have deep radio and X-ray 
observations which may identify starbursts and AGN. Mag- 
nitudes in this paper are given in the Vega system and can 
be converted to the AB system using Rab = Rvega + 0.19 
and Kab = Kyega + 1.87. Quantities dependent on Hub- 
ble's Constant are given in terms of /2.50 = Ho/50 or hioo = 
Ho/100 km s _1 Mpc _1 . 

This paper is organized as follows: Section 2 describes 
the dataset and its reduction and analysis, and Section 3 the 
selection of the EROs and their distribution on the sky. In 
Sections 4 and 5 we investigate, in turn, the number counts 
and angular correlation function, ui(6), of the K < 21.0 
EROs and other galaxies, and compare with models. In 
Section 6 we fit radial profiles to a subsample of brighter, 
K < 19.5, EROs to investigate their morphologies and an- 
gular sizes, and in Section 7 describe the radio properties. 
8 discusses these findings and with other recent studies of 
EROs , and their possible interpretation. 



2 OBSERVATIONS 

2.1 Data: K and R Imaging 

Our K (A = 2.2/im) and R band imaging covers an area 
within field N2 of the European Large Area ISO Survey 
(ELAIS, Oliver et al. 2001), centered at R.A. 16 h 36 m 30 a 
Dec +41:04:30. The first part of the if-band data consists 
of a mosaic of 16 contiguous fields observed using the UKIRT 
Fast-Track Imager (UFTI), on Mauna Kea, Hawaii, between 
1 January and 31 May 2000. Most frames received a total 
of 8000 seconds exposure time. The UFTI camera contains 
a 1024 x 1024 pixel HgTeCd array which, with a pixel size 
of 0.091 arcsec, covers 1.55 x 1.55 arcmin. 

The second part consists of three fields observed with 
the Isaac Newton Group Red Imaging Device (Ingrid) on the 
William Herschel Telescope (WHT), La Palma. The Ingrid 
fields lie on the edges of the UFTI mosaic, overlapping it 
slightly. Ingrid is fitted with a 1024 x 1024 pixel near-IR 
detector with a pixel size of 0.238 arcsec, covering 4.06 x 4.06 
arcmin. Our Ingrid data are not quite as deep as the UFTI 
mosaic, due to a shorter exposure time (5540 sec) and lower 
instrumental sensitivity. 

The UFTI and Ingrid observations lie within the area 
of an 7?-band image obtained with the Prime focus cam- 
era on the WHT (May 1999). This instrument contains 2 
mosaiced 2048 x 4096 pixel EEV chips, with a pixel size of 
0.238 arcsec, covering 16 x 16 arcmin. The i?-band data, 
obtained in May 1999 as part of a study of faint Chan- 
dra X-ray sources (Gonzalez-Solares et al. 2002), consists of 
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13 spatially dithered 600 sec exposures, combined to give 
an octagonal 16.4 x 16.9 arcmin frame with detection limit 
R ~ 26. This has accurate astrometry (to < 0.5 arcsec) 
based on the positions of radio sources. 

2.2 Radio/X-ray /sub-mm Data 

The ELAIS N2 field has been surveyed in a number of other 
passbands. Radio (VLA) observations (Ivison et al. 2002) 
reach a 3a detection limit of F(1A GHz) = 30/iJy (with 
beam size 1.40 x 1.46 arcsec), and detect > 100 sources. 
The available X-ray data consist of a 75ks Chandra ob- 
servation made on 2 August 2000, covering a 16.9 x 16.9 
arcmin field, in which 91 sources were detected to J^(0.5- 
8.0keV) ~ 6 x 10 -16 ergs s^cm -2 (Manners et al. 2002; 
Almaini ct al. 2002), most of which now have optical iden- 
tifications (Gonzalez-Solares et al. 2002). Sub-mm observa- 
tions with SCUBA, detected 17 sources above a 3.5cr limit 
F(850^m) ~ 8 mjy (Scott et al. 2002). 

2.3 Data Reduction 

The UFTI and Ingrid data were reduced and analysed us- 
ing I RAF. The UFTI dataset was made up of 16 pointings, 
each observed for 9 spatially dithered exposures of 800 sec 
(two pointings had slightly different exposures of 700 sec 
and 1100 sec and were renormalized to 800 sec). Our in- 
tention was to combine all 9 x 16 = 144 exposures into a 
single mosaiced image, and to this end, the spacing between 
the grid of pointings had been matched to the field-of-view 
and dither pattern so that the resulting mosaic would have 
an approximately uniform 7200 sec coverage (except at the 
edges) . 

The sky background was subtracted from all UFTI data 
(using IRAF 'sky'). Astrometry was derived for each point- 
ing, by reference to detections on the WHT R image. How- 
ever, even with these transforms, the raw exposures could 
not be mosaiced with acceptable accuracy, indicating slight 
rotations or other distortions between them. This problem 
was remedied by choosing one pointing as a reference and 
then using IRAF 'wcsmap' to fit a general polynomial trans- 
form between its co-ordinate system and that of the other 15. 
These 15 frames were then rebinned, using IRAF 'geotran', 
into the co-ordinate system of the reference. The rebinned 
data could then be simply combined, using IRAF 'combine' 
(with 'offsets=wcs'), into a 16- frame UFTI mosaic (taking 
care to exclude any 'bad' regions from each exposure in 
the summation). The final step was to trim the noisy edges 
(which received less than the full exposure time) from the 
mosaic, to leave a total usable area of 38.7 arcmin . 

The three fields of Ingrid data overlap slightly with the 
UFTI mosaic, but not with each other, and were reduced 
separately. The data for each field consisted of 90 spatially 
dithered exposures of 61.56 seconds (each in turn made up 
of 6 exposures which were co-added during observing). All 
exposures were debiased, and flat-fielding was performed for 
blocks of 9 consecutive exposures, by dividing them by their 
median image (derived using 'combine' with sigclip rejec- 
tion, and normalized to a mean of unity). 

The spatial dithers, which followed the same sequence 
for each block of 9 exposures, were measured using 'xregis- 
ter'. All 90 exposures of each of the three fields could then 



be added using 'combine'. Magnitude calibration was de- 
termined from short exposures of UKIRT standard stars, 
interspersed with the observations. We then matched ~ 10- 
15 stars on each Ingrid field (detected as described below) to 
their known RA and Dec on th WHT R image and fitted as- 
trometric transforms using IRAF 'pltsol' (with rms residuals 
0.11-0.25 arcsec). 

2.4 Source Detection 

Sources were detected on the K and R images using SEx- 
tractor (Bertin and Arnauts 1996). Throughout, magnitudes 
given are the 'total' magnitudes derived by SExtractor by 
fitting elliptical apertures to each individual detection. 

For the UFTI mosaic, our detection criterion was that 
a source must exceed a 1.5a s k y threshold (21.36 K mag 
arcsec" 2 ) in 16 contiguous pixels (0.13 arcsec 2 ), and we also 
used a detection filter of Gaussian FWHM 4 pixels. Using 
a plot of detection magnitude against FWHM, stars could 
be separated from galaxies to K = 16.5 and many spuri- 
ous noise detections could be rejected. The mean FWHM 
of unsaturated stars was only 0.67 arcsec, indicating very 
good seeing, although there is significant variation in the 
PSF between the different frames of the mosaic. 

For the WHT _R-band image, our detection threshold 
was 1.5a B k y (26.66 R mag arcsec -2 ) in a minimum area 4 
pixels (0.23 arcsec 2 ), with a detection filter of 2.0 pixels 
Gaussian FWHM. The total usable area is 212 arcmin 2 . The 
mean stellar FWHM is 0.76 arcsec, and from the number 
count of detected galaxies, detection appears to be complete 
to R ~ 25.5, with moderate (25-30 per cent) incompleteness 
at 25.5 < R < 26 and a turnover at R > 26. The Ingrid 
data (total area 49.5 arcmin 2 ) have a similar pixel size and 
again we used the detection criterion of 1.5a s k y (20.60-20.65 
R mag arcsec -2 ) in > 4 pixels, with a Gaussian 2.0 pixel 
FWHM filter. Stars (45) were separated from galaxies to 
K = 16.5. The mean stellar FWHM of 0.68 arcsec again 
indicated good seeing. The galaxy counts from the Ingrid 
data are very similar to the UFTI counts to K = 19.5 but 
turn over at K > 20 rather than at K > 21, hence are ~ 1 
mag less deep. 

The next step was to merge the Ingrid and UFTI cat- 
alogs, by flagging for exclusion the Ingrid detections within 
the area already covered by the UFTI mosaic. The Ingrid 
and UFTI magnitudes of the overlap- region objects were 
generally consistent within ~ 0.1 mag. Excluding overlaps 
reduced the total catalog area to 83.6 arcmin 2 . Hereafter we 
make use of the combined UFTI+Ingrid catalog to K — 20 
and the UFTI mosaic only at 20 < K < 21. 



3 IDENTIFYING THE R - K > 5 GALAXIES 

3.1 The Significance of R K > 5 

Galaxies with observer-frame R — K > 5.0 (R — K > 3.31 
in the AB system) are either old, passive (zero or near-zero 
SFR) galaxies at z > 0.9, or are very dust-reddened. Figure 
1 shows R— K against redshift for three models, computed 
using Pegase2 (Fioc and Rocca-Volmerange 1997) with the 
time-redshift relation for Q m = 0.3, S7a = 0.7 and hioo = 
0.55. 
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Figure 1. Observed R—K colour against redshift, for three mod- 
els representing an evolving elliptical galaxy (pERO), a heavily 
reddened starburst (dsfERO), and a normal spiral (which is never 
an ERO). 

(i) An evolving model representing an elliptical galaxy, 
which forms stars in a 1 Gyr dust-reddened burst at z > 3.4. 
The SFR and dust then fall to zero, and passive evolution 
occurs. The galaxy would have R — K > 5.0 at all z > 0.93 
and hence is a 'pERO'. 

(i) A pure, non-evolving starburst, with very strong 
dust reddening of E(B - V) = 1.0 mag. Such 'dsfEROs' 
are reddest at 1 < z < 1.6, but if as dusty as this model can 
be EROs at all 0.5 < z < 3. Starburst galaxies with an old 
stellar component, or post-starbursts, could be EROs with 
somewhat less reddening. 

(iii) For comparison, a model representing an evolv- 
ing spiral (Sbc) with near-solar metallicity and a 'normal' 
amount of dust. 'Normal' star-forming spirals like this do 
not have R — K > 5 at any redshift. 

3.2 ERO Selection 

The /^-detected galaxies were matched in RA and Dec with 
detections on the WHT i?-band image, within a maximum 
offset radius of 1.5 arcsec (two small areas of Ingrid data 
were not covered in R, slightly reducing the K + R overlap 
region to 81.5 arcmin 2 ). The colour of each galaxy was taken 
to be simply the difference of the R and K total magnitudes. 
To test this, the 'total mag' colours were compared with 
R — K measured in fixed circular (2.5 arcsec) apertures, for 
a subsample of EROs. These colours agreed within < 0.15 
mag in most cases, with no significant systematic difference. 

Any K < 21 detection with either R — K > 5 or no 
R counterpart was nagged as a candidate ERO. The K de- 
tections in the latter category are either (i) spurious, (ii) 
real galaxies where the R counterpart is obscured (e.g. by a 
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Figure 2. R—K colour against K for UFTI and Ingrid detections 
to K = 21, with star-galaxy separation to K = 16.5. Galaxies 
undetected in R, and hence included as EROs, are placed at R = 
30, but may lie anywhere above the R = 26 locus. 

diffraction spike) or merged with a brighter detection, (iii) 
real galaxies too faint in R to be detected on our image, 
which can reasonably be classified as EROs. All were care- 
fully examined by eye, on the K image and at the corre- 
sponding position on the R image. In this way, a number of 
spurious detections were identified, and thereafter excluded. 

Our final catalog of verified EROs consists of 158 ob- 
jects (including 49 undetected in R), of which 99 are on the 
UFTI mosaic. Figure 2 shows R — K against K magnitude 
for the detected galaxies - EROs appear only at faint mag- 
nitudes, mostly K > 18. One much brighter (K = 16.2) 
object has R — K = 5.72, but appears stellar and is almost 
certainly a red Galactic star, and is not included here as an 
ERO. 

3.3 Spatial Distribution 

Figure 3 shows the distribution of EROs. ERO maps are 
useful for identifying galaxy clusters at z > 1, where the 
early-type members will show up as an arcmin-scale over- 
density (e.g. Stanford et al. 1997). In our data, there may 
be a group of EROs centered on UFTI detection number 
608, a bright ERO (K = 17.78) of particular interest in that 
it is a radio and X-ray source (Section 7). Within a 45 arcsec 
radius (shown on Figure 3) of its position there are 8 EROs 
with K < 20, compared to 2.4 expected for a random dis- 
tribution. These EROs are discussed further in Section 6.3. 



4 NUMBER COUNTS OF GALAXIES AND 
EROS 
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Figure 3. Distribution on the sky of the K < 20 EROs (R-K > 
5.0) on the UFTI mosaic (central area) and the 3 Ingrid fields, 
with the 20 < K < 21 EROs on the UFTI field only. The plotted 
circle shows a 45 arcsec radius around UFTI detection 608. 

4.1 Observations 

Figure 4 shows if-band differential number counts for all 
galaxies, from the combined UFTI+Ingrid data at K < 20 
and UFTI data only at 20 < K < 21. Also shown are counts 
from the deeper (to K ~ 23) Keck and ESO-VLT surveys 
of Moustakas et al. (1997) and Saracco et al. (2001). Our 
counts are reasonably consistent with these, although lower 
at K > 20, suggesting some incompleteness - the ratio of 
our galaxy count to that of Saracco et al. (2001) falls from 
0.968 at 19.5 < K < 20.0 to 0.725 at 20.0 < K < 20.5 and 
0.611 at 20.5 < K < 21.0. 

Figure 5 and Table 1 show number counts for EROs 
only. The fraction of EROs in the UFTI+Ingrid data is 
112/812 = 13.8 per cent at K < 20 and 158/1076 = 14.7 
per cent at K < 21.0. The 'incompleteness corrected' count 
is derived by dividing our observed ERO count by the ra- 
tio of our count for all galaxies to that of Saracco et al. 
(2001) at the same K limit. Also plotted are the R~K > 5 
ERO counts from Daddi et al. (2000), and the Saracco et al. 
(2001) ERO counts - note that the latter are selected with 
a different criterion of J — K > 1.9, which may include a 
greater number of 'dsfEROs'. Our ERO counts agree well 
with these two surveys. 

4.2 Models: PLE, NE, Merging and 'M-DE' 

These counts are compared with simple galaxy evolution 
models, with and without the effects of merging. For the 
models considered in this paper, we adopt cosmological pa- 
rameters of Qm = 0.3, S!a = 0.7 (using the analytic lu- 
minosity distance formula of Pen 1999) and the obserav- 
tional (from 2MASS) K-b&nd galaxy luminosity functions of 



Figure 4. ii"-band differential galaxy number counts, for galax- 
ies of all colours, from our data (Ingrid+UFTI to K = 20, UFTI 
only at 20 < K < 21), and the deeper surveys of Moustakas et 
al. (1997) and Saracco et al. (2001), compared with three mod- 
els; non-evolving (long-dashed), PLE (dotted), and merging with 
R<j> = Rm* = 0.3 (short-dashed). 



Table 1. Observed number counts (number N g and surface den- 
sity, p) of R- K > 5.0 EROs on the UFTI+Ingrid field, with y/N 
errors, and (at K > 19.5) p with an estimated incompleteness 
correction. 



K 


N 3 


p (Obs.) 


p ('Corrected') 


interval 




deg _2 mag _1 


deg _2 mag _1 


17.0-17.5 


1 


88.4+ 88.4 




17.5-18.0 


3 


256.0+ 153.0 




18.0-18.5 


8 


706.8 ± 249.9 




18.5-19.0 


23 


2032.2 + 423.7 




19.0-19.5 


38 


3357.5 + 544.7 




19.5-20.0 


39 


3445.9 + 551.8 


3560. 


20.0-20.5 


20 


3827.2 ± 855.8 


5279. 


20.5-21.0 


26 


4975.4 ± 975.8 


8143. 



Kochanek et al. (2001). luminosity evolution is modelled us- 
ing Pegase2 (Fioc and Rocca-Volmerange 1997) with a range 
of star- formation histories, and (as in Roche et al. 2002), an 
initial mass function with the Salpeter slope x = 2.35 (where 
|^ oc ]VT X ) at 0.7 < M < 120 M Q , flattening to x = 1.3 at 
0.1 < M < 0.7 M Q . 

Galaxies in the early-type luminosity function (repre- 
sented as a Schechter function with = —25.04, a — 
-0.92 (f>* = 0.0005625 Mpc" 3 for h 50 = 1) form all their 
stars in an initial dust-reddened (E(B — V) = 0.65 mag ) 
burst, beginning 16 Gyr ago (z ~ 6), with a range of du- 
rations represented as either 1 Gyr (for half) or 2 Gyr. Af- 
ter the burst both dust and SFR fall to zero and evolution 
is thereafter passive. Galaxies in the late- type luminosity 
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gives f- = -t H (l + 2) 2 ■ Substituting, 
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Figure 5. _R"-band differential galaxy number counts for R — 
K > 5.0 EROs from our UFTI+Ingrid data and Daddi et al. 
(2000), and for J - K > 1.9 EROs in Saracco et al. (2001), with 
minimal ,/N errors, compared with R — K > 5.0 galaxy counts 
predicted by PLE (dotted), non-evolving (long-dashed), merging 
with R^ = Rm = 0.3 (short-dashed), and 'M-DE' merging with 
Rm = 0.3 R,p = —0.46 (dot-dashed) models. 



function (M* K = -24.49, a = -0.87 4>* = 0.0012625 MpcT 3 
for /Y50 = 1) are represented with a range of continuous 
star-formation models with different timescales and moder- 
ate dust (in the amounts calculated by Pegase). The only 
R — K > 5 galaxies in this model are the E/SOs at 2 > 0.93. 

In the pure luminosity evolution (PLE) model the 
galaxies evolve only in L* , with no change in <f>* or a. As our 
version of this model gives relatively strong L* evolution, we 
also show a non-evolving (NE) model, with the same lumi- 
nosity function at 2 = 0. In our merging model, the evolution 
of luminosity per unit mass in the same as the PLE model, 
but superimposed on this is an increase in comoving number 
density (0*) and an associated decrease in the characteristic 
galaxy mass (m*) with redshift, cancelling out some of the 
evolution in galaxy luminosity (L*). 

The merger rate is parameterized in terms of the effect 
of merging on galaxy mass, A ^ ^ = R m (z) per Hubble time 
(tn)- Observationally, Patton et al. (2001) estimated from 
CNOC2 data that R m (z) ~ 0.3(1 + z) 2 ' 3 to 2 = 0.55. If 
we assume this evolution to continue to z = 1, with R m (z) 
constant at z > 1, then m*(z) will evolve as: 



dm* „ , . 1 , .2.3 
= m z R m (0) t H (1 + 2) 



at 2 < 1 



dm* * , . _i_2.3 , ^ 
—jj- = m z R m (0) t H 2 at 2 > 1 

approximating the lookback time t n ow~t ~ tu — t£r(l+2)~ 1 , 



^f- = -m* z R m (0) t^il + zf 3 at2<l 
az 

^- = -m* z R m (0) ^ 1 2 2 ' 3 (l + 2)- 2 at2>l 



Integrating, 



[ d ^-^)f\^ z rdz 

In m*(z) = In m*(0) - R m (0)((l + 2) 1 ' 3 - 1)/1.3 

m*(z) = m*(0 exp [-i? m (0)((l + z) 1 " 3 - 1)/1.3] 2 < 1 
and, 



(l + zy 2 dz 



In m*(z) = In m*(l) - 2 2 ' 3 fl m (0)(i - (1 + z)' 1 ) 

m*(z) =m*(l) exp [-2 1 ' 3 i? m (0)(l - 2(1 + z)' 1 )] at 2 > 1 

with R m (0) = 0.3 this 

gives m*(l,2,3) = (0.714, 0.558, 0.493)m* (0). If the corre- 
sponding opposite evolution in comoving number density, 
R4,, occurs at the same rate (R4, = R m ), galaxies at 2 ~ 3 
are half as bright but twice as numerous as they would be 
in a PLE model. 

The counts of all galaxies (Figure 4) exceed the NE 
model, and are more consistent with PLE or merging. The 
observations do not firmly distinguish between these two, 
as the model counts are relatively insensitive to merging. 
Beyond the limits of our survey there is some indication 
that the count exceeds all three models, which may indicate 
the luminosity function has a steeper faint-end slope than 
derived by Kochanek et al. (2001). 

On Figure 5 the same models are used to predict the 
counts of R-K > 5.0 EROs. Although only the 'pERO' type 
is included, our PLE model overpredicts the ERO counts at 
all magnitudes, by factors ~ 3. The inclusion of R^ — R m — 
0.3 merging improves the fit at the bright end, but increases 
the overprediction at K > 20. However, the non-evolving 
model does fit the ERO counts. 

Firth et al. (2002) similarly found number counts of 
EROs to be well-fitted by a non-evolving model to H = 
20 (K ~ 21), and also claim that this model fits their 
photometrically-estimated redshift distribution N(z). Our 
PLE model, with its relatively top-heavy IMF and ejection 
of the dust at the end of the starburst, gives strong L* evolu- 
tion for E/SOs, and in combination with the high formation 
redshift, predicts a high ERO count. 'Milder' PLE models 
may be somewhat closer to the observed ERO counts, but 
to be consistent with them, we and Firth et al. (2002) had 
to reduce L* evolution to near-zero. Cimatti (2000b) did fit 
the ERO counts with a PLE model, but only for an ellipti- 
cal formation redshift 2 < 2.2, which is unrealistically low 
in view of the existence of Lyman break objects, SCUBA 
sources and radio galaxies at 2 > 3. Hence we conclude that 
the number counts of EROs are in general inconsistent with 
PLE models. 
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Figure 6. Features of the 'M-DE' model; (above) evolution of 
characteristic galaxy mass m* and comoving number density <p* ; 
(below) evolution of the characteristic luminosity L* , in the rest- 
frame K and R bands, and of the luminosity per unit mass, for an 
elliptical with a 1 Gyr initial starburst (star-formation at 3.4 < 
z < 6). 



With the aim of finding a physically realistic evolv- 
ing model for ERO counts, we consider the merging model 
again, and retain R m — 0.3 while varying to best-fit our 
incompleteness-corrected ERO count. The minimum x 2 oc- 
curs for R,), ~ —0.46 (with la error bar ±0.10). Figure 5 
shows the counts for this model, hereafter 'M-DE' (merging 
with negative density evolution). They are very similar to 
the NE model counts, to K ~ 20, and in good agreement 
with observations. 

Figure 6 illustrates the features of the M-DE model. 
Both mass and comoving number density decrease with red- 
shift. Luminosity per unit mass (shown for restframe R 
and K) increases with redshift due to passive evolution (at 
z > 3.4 the galaxy is starbursting and most of the red light 
is hidden by dust). Luminosity (L*) increases with z at a 
slower rate as it is partially cancelled by the decrease in 
mass. Physically, this model could be represented as 60 per- 
cent of the present-day comoving number density of E/S0 
galaxies forming at z > 1, and 42 percent at z > 2, with 
the remainder originating at lower redshifts, e.g. from merg- 
ers of disk galaxies in which the remaining gas content is 
exhausted in a merger-triggered starburst. Some spiral-to- 
elliptical transformation could also occur through the cumu- 
lative effect of non-merging interactions, or gas-stripping in 
clusters. 

In addition to rejecting the PLE model, the ERO counts 
do not support the other extreme case of 'pure' hierarchical 
merging. Cimatti (2002b) compares the observed count of 
passive EROs at K < 19.2 with the predictions of two such 
models from Firth et al. (2002) and Smith et al. (2002), both 
derived from ACDM cosmologies. These models predict, re- 



spectively, 0.19 and 0.04 EROs arcmin -2 at K < 19.2. 
Our observed (UFTI+Ingrid) count of EROs to this limit 
is 0.60 ± 0.09 arcmin~ 2 , and hence a mixture of formation 
histories, as represented by the M-DE model, is preferred. 



5 CLUSTERING OF EROS 

5.1 Calculating the Angular Correlation Function 

The clustering properties of the EROs will provide clues to 
their nature and history. We investigate the clustering by 
calculating the angular correlation function, to(9), for both 
R — K > 5 galaxies and full AT-limited samples, by the meth- 
ods described below. To K — 20 we use the combined UFTI 
and Ingrid catalogs, but for fainter limits the analysis is lim- 
ited to galaxies on the UFTI mosaic. 

For N g galaxies brighter than a chosen magnitude limit, 
there will be ^N g (N g — 1) possible galaxy-galaxy pairs. 
These are counted in bins of separation of width A (log 9) = 
0.2, giving a function N gg (9i). A large number of random 
points (N r = 50000 here) is scattered over the same area 
as covered by the real galaxies, and the separations of the 
N g N r galaxy-random pairs, taking the real galaxies as the 
centres, are counted in bins to give N gr (9i). The separations 
of the ^N r (N r — 1) random-random pairs are also counted 
to give N rr (9i). 

If DD = N gg (9i), and DR and RR are the galaxy- 
random and random-random counts normalized to have the 
same summation over 9, 



DR ■- 



RR = 



2N r 



-N gr (9i) 



-JVrr(0i) 



N R (N r - 1) 

then, using the Landy and Szalay (1993) estimator, 
DD - 2DR + RR 



u>{9i) 



RR 



Errors were estimated by dividing the data area into 20 sub- 
areas and then recalculating co(9) for 20 subsamples, each 
time excluding both the real galaxies and the random points 
from a different sub-area and using the remaining 19. The 
scatter between the lo(9) of these 20 subsamples is then mul- 
tiplied by ^ 19 20 9 — 4-25 to give the error bars for the full 
dataset uj(9). 

This estimate will be negatively offset from the true 
to(9) due to the restricted area of observation (the 'integral 
constraint'). If the real u>(9) is of the form A LJ 9~ S , where 
is an amplitude, the estimate corresponds to A LJ (9^ S — C). 
The negative offset AC can be estimated by doubly inte- 
grating an assumed true ui(9) over the field area SI, 



AC 



= *JJ 



uj(9)dn 1 dn 2 



Using the random-random correlation, this can be done nu- 
merically - 



C = 



Y,Nrr{9)9~ 



T.N rr {9) 

Assuming 8 — 0.8, this gives C = 9.226 for the combined 
UFTI+Ingrid area, and 13.46 for the UFTI mosaic only. 
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Figure 7. Angular correlation functions, u)(9) as calculated for 
R — K > 5.0 EROs, and for all galaxies, on the combined 
UFTI+Ingrid field to limits of (a) K < 19.5 and (b) K < 20.0, 
with best-fit functions of the form 'A u (6-°- 8 - C)' (dotted). 



The amplitude A w is then obtained by least-squares 
fitting A(0-°- s - C) to the observed w(0), over the range 
1.26 < < 502 arcsec, weighting each point using the error 
bars estimated as above. To estimate an error on A w , the 
same function is then fitted to the uj(6) of each of the 20 sub- 
samples, and multiplying the scatter between the subsample 
A w by, again, 4.25. 



5.2 w(0) results 

For the EROs we detect positive clustering at the > 2a 
level at magnitude limits K = 19.5-20.0. To K = 20.25 the 
smaller sample provided by the only the UFTI data gives 
only 1.5a evidence of clustering. At K = 20.5-21.0, there is 
no detection even at lcr, but we do not consider these limits 
here as w(0) is likely to be affected by the incompleteness. 
For the full if-limited sample the detection of clustering 
never reaches 2a, and the upper limits indicate the A u is 
about an order of magnitude lower than for EROs at the 
same K limits. 

Figure 7 shows the observed w(0), with the fitted func- 
tions, at two Jf-band limits, and Table 2 gives the best-fit 
Au with ±1it errors. Figures 8 and 9 show our A^ against K 
limit, together with results from previous clustering analyses 
and models (described below). In general our results are con- 
sistent with all the published A^ , support the strong cluster- 
ing of EROs reported by Daddi et al. (2000) and Firth et al. 
(2001), and suggest that this continues to at least K ~ 20. 



Table 2. Galaxy u)(9) amplitudes A u (in units of 10 4 at one 
degree) of full A"-limitcd galaxy samples and of EROs (R — K > 
5.0), as estimated from UFTI+Ingrid data (UFTI only at K > 
20). N g is the number of galaxies in each sample. 



K mag 


All 


galaxies 


R - 


K > 5.0 EROs 


limit 


N g 


Au 


N g 


Au 


19.50 


615 


9.30 ± 5.26 


73 


151.2 ± 59.2 


19.75 


714 


5.32 ±4.48 


93 


86.63 ±35.26 


20.00 


813 


5.76 ±3.65 


112 


76.61 ± 33.25 


20.25 


418 


8.24 ± 4.93 


63 


36.02 ± 24.29 



"T" 



"T" 



B Baugh et al. (1996) NEP field 
□ Baugh et al. (1996) NGP field 
A Roche et al. (1998) 
Q Roche et al. (1999) 
Kuemmel and Wagner (2000) 
A Daddi et al. (2000) 
ifMcCracken et al. (2000) 
• This work 



merg. e = — 1,2 r =5.85 



\NE e=0 r =5.85 
merg>e-^0^r =5.85 

PLE E = - r 



J I L 

20 

K Magnitude Limit 



_J 

22 



24 



Figure 8. The scaling of u>(8) amplitude with K magnitude limit 
for full A"-limitcd samples of galaxies, as derived from our data, 
shown with previously published results and four models (all with 



Hm = 0.3 Qa = 0.7 and ro 
e = (dotted), merging (with e 



5.85 ftToo Mpc) - PLE with 



(short-dash), NE with e = 
(long-dash) and M-DE with comoving (e = —1.2) clustering, (dot- 
dashed). 



5.3 Modelling and Interpretation of uj{6) 

The observed A^ of any sample of galaxies depends on their 
redshift distribution N(z), and the cosmological geometry 
(through the angular size distance), as well as on the in- 
trinsic clustering of the galaxies, as parameterized by their 
two-point correlation function in three dimensions, £(r, z). 
Hence to interpret A u , we begin by modelling N(z). 

Figure 10 shows the N(z) of EROs, at K < 19.5 and 
K < 21.0, as given by our NE, PLE and M-DE models. 
PLE gives the highest numbers of EROs and most extended 
N(z), while the M-DE model N(z) is much lower in nor- 
malization and a little less extended in redshift, especially 
at the brighter limit. The NE model, as previously noted, 
gives similar numbers of EROs to M-DE, but the N(z) is 
even less extended and peaks at a lower redshift. 

The intrinsic (three-dimensional) clustering of galaxies 
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Figure 9. As Figure 8 for EROs (R-K > 5.0 galaxies) only. The 
//-band limits of Firth et al. (2002) are converted to K assuming 
H — K = 1 (approximately valid for 0.8 < z < 2.6 ellipticals). 
The plot shows four models fitted to the observed ui(6) scaling; 
comovmg M-DE with r = 12.8 h~^ Q Mpc (dot-dash), stable M- 
DE with ro = 21.4 ^ u) o Mpc (long-dash), comoving NE with ro = 
10.3 Mpc (short-dash) and stable NE with ro = 16.4 h^Q 

Mpc. Also plotted, for comparison, is the M-DE stable model with 
ro = 5.85 h^Q Mpc, showing the order-of-magnitude difference 
between ERO cluatcring and that of 'average' faint galaxies. 



and its evolution with redshift is represented here by the 
power-law form 

t{r,z) = (r/r a )- 1 (l + z)-^ ) 

where ro is the local correlation radius, 7 ~ 1.8 (obser- 
vationally) and e represents the evolution with redshift, 
e = being clustering stable in proper co-ordinates, and 
e — —1.2 comoving clustering. This produces a projected 
(two-dimensional) clustering, u>(9) = A^f? - ' 7-1 ', where the 
amplitude is given by Limber's formula (see e.g. Efs- 
tathiou et al. 1991), 

-(3+e) 



Jo 



, dx(z) 
dz 



■[(N(zY]dz/[ 



I 



N(z)dzY 
3.679 for 7 



xi~ 1 (z) 

where x(z) is the proper distance and C 7 
1.8. 

Figure 8 shows four models for the u>(9) scaling of 
all galaxies with A"-band limit, all with the normalization 
ro = 5.85 /ij~ g Mpc, as estimated by Cabanac, de Lappar- 
ent and Hickson (2000) from an /-band survey. The NE 
model predicts a higher A u than the PLE and merging 
(Rtf, — R m * = 0.3) models, due to the less extended N(z), 
and a negative e increases A w . 

The observed A w are consistent with or slightly above 
the merging and PLE models with e ~ 0. The NE and es- 
pecially the comoving model overpredict the lj(6) to some 
extent. 



Figure 10. The redshift distribution N(z) of R—K > 5.0 galaxies 
(pEROs) predicted by PLE (dotted), NE (long-dashed) and M- 
DE (dot-dash) models, at limits K < 19.5 and K < 21.0. 



Figure 9 shows some models which fit the A u of EROs. 
Giant ellipticals are found in local surveys to be more 
strongly clustered than 'average' galaxies; e.g. Guzzo et al. 
(1997) estimate r = 8.35 ± 0.75 h^ Mpc for E/SOs with 
Mb < —21. If the EROs, or most of them, are indeed the 
progenitors of massive ellipticals, they would be expected 
to show similarly strong clustering. To estimate the intrin- 
sic clustering (ro) of EROs, we very the ro to minimize the 
X 2 for the 7 plotted data points (from all three surveys) at 
18.8 < K < 20.0, for our M-DE and NE models, and for 
both e = and e = -1.2. 

For comoving clustering, we best-fit r = 12.8T1.5 h^ 



Mpc for the M-DE model and r 
NE, and for stable clustering, ro : 



: 10.3 ±1.2 Mpc for 
21.4 ±2.6 him M P C for 



the M-DE model and r = 16.4 ± 2.0 h~ l Mpc for NE. 
We note that Firth et al. (2002) estimate very similar ro 
for H < 20.0 EROs and a photometrically estimated N(z). 
Our four fitted models are very similar to K = 21, and all 
are consistent within la of the observations - much deeper 
data will be needed to constrain the evolution within the 
-1.2 < e < range. 

In summary, for the probable range of E/S0 luminosity 
evolution, we estimate the intrinsic clustering of bright, K < 
20, EROs, in the form of a comoving correlation radius, as 
ro — 10-13 fej'oo. This implies that the intrinsic clustering of 
at least the brighter EROs is indeed stronger than that of 
'average' galaxies, and even exceeds the ro ~ 8ftJ" o Mpc of 
present-day giant ellipticals, although it may not quite equal 
the ro — 21/i 1 ~ 1 Mpc of local Abell clusters (Abadi, Lambas 
and Muriel 1998). This is discussed further in Section 8.2.2. 
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Table 3. Co-ordinates (equinox 2000.0) of the 32 EROs in our bright sample (K < 19.5 on UFTI image), K and R magnitudes, best-fit 
r^i (in arcsec, with error) ellipticity 1 — ^ (0 = round), best-fit profile type (d = disk b = bulge p = point-source), significance A(x 2 ), 
in units t(x 2 ) of favouring disk or bulge, point-source fraction f p , and other properties of morphology (int = in interacting pair, dn = 
double nucleus, asn = asymmetric - offcentre nucleus, aso = asymmetric outer regions, ir = irregular/flocculcnt) ) or emission (X = 
Chandra x-ray detection, Ra = VLA radio detection). 



Number 


ti.A 


Dec. 


K 


R 


rhl 


ell 


Type 


sig 


fp 


other 


r91 


16:36:54.32 


40:56:57.43 


19.084 


25.655 


0.35 ±0.10 


0.188 


b 


2.1 







rl96 


16:36:47.92 


40:57:52.13 


19.034 


25.859 


0.18 ±0.04 


0.128 


b 


1.03 







r256 


16:36:48.43 


40:58:08.46 


18.881 


25.038 


0.38 ±0.15 


0.061 


b + p 


2.1 


0.375 


Ra 


r270 


16:36:45.13 


40:58:15.69 


18.843 


24.625 


0.28 ±0.04 


0.229 


b 


4.0 







r280 


16:36:47.87 


40:58:22.07 


19.325 


24.376 


0.20 ±0.03 


0.171 


b 


1.05 







r332 


16:36:50.31 


40:58:42.15 


19.226 


24.759 


0.24 ±0.03 


0.118 


d 


0.47 





asn 


r407 


16:36:46.62 


40:59:16.15 


19.472 


> 26 


0.78 ±0.28 


0.106 


d + p 


3.3 


0.6 


int 


r492 


16:36:53.51 


40:59:31.85 


19.309 


25.951 


0.19 ±0.07 


0.019 


b + p 


0.38 


0.5 




r506 


16:36:42.77 


40:59:34.36 


19.410 


24.997 


0.31 ± 0.09 


0.090 


b 


0.05 





irr 


r518 


16:36:53.47 


40:59:37.42 


19.132 


> 26 


0.59 ±0.18 


0.230 


b 


3.1 





aso Ra 


r552 


16:36:56.25 


40:59:47.74 


19.194 


25.288 


0.32 ±0.04 


0.039 


b 


0.05 





Ra 


r561 


16:36:53.89 


40:59:55.92 


18.019 


23.246 





0.043 


V 


- 


1.0 




r581 


16:37:00.60 


41:00:06.62 


19.356 


26.367 


0.34 ±0.03 


0.248 


d 


0.10 







r585 


16:36:46.19 


41:00:05.93 


19.131 


26.950 


0.40 ±0.10 


0.248 


b 


1.67 





asn 


r594 


16:36:39.43 


41:00:07.48 


19.377 


26.338 


0.32 ±0.10 


0.210 


b + p 


0.81 


0.25 


aso 


r599 


16:36:45.46 


41:00:10.81 


19.463 


25.957 


0.42 ±0.14 


0.062 


b 


0.08 





asn/ aso 


r608 


16:36:57.76 


41:00:21.71 


17.780 


23.176 


0.50 ±0.06 


0.152 


b + p 


26.1 


0.09 


X Ra 


r622 


16:36:45.80 


41:00:23.72 


18.501 


24.802 


0.66 ±0.05 


0.325 


d 


14.8 





dn 


r626 


16:36:58.36 


41:00:25.34 


19.338 


25.311 


0.22 ±0.07 


0.056 


b 


0.67 









1U.OU.OO.O-L 


/ii -nn-98 1 1 

11.UU.iO.il 


1 fi 798 

lO. / iO 


9zL fi^R 




n 1 rizi 

U. 1U1 


u 


9 R 
z.o 


n 
u 




r642 


16:36:56.99 


41:00:31.76 


18.610 


24.819 


0.32 ±0.04 


0.105 


b 


7.9 





Ra 


r650 


16:36:38.92 


41:00:33.62 


18.984 


24.466 


0.59 ±0.07 


0.438 


d 


16.4 





irr 


r660 


16:36:41.05 


41:00:38.86 


19.248 


25.062 


1.42 ±0.30 


0.582 


d 


12.9 





irr Ra 


r675 


16:36:56.29 


41:00:46.77 


18.695 


23.960 


0.61 ± 0.05 


0.429 


d 


> 20 







r852 


16:36:49.75 


41:01:32.89 


19.222 


25.575 


0.28 ±0.14 


0.173 


b + p 


1.79 


0.10 




r952 


16:36:45.31 


41:02:12.34 


19.111 


> 26 


0.57 ± 0.08 


0.351 


d + p 


> 2 


0.14 


Ra 


r955 


16:37:02.38 


41:02:19.25 


17.471 


22.687 


0.38 ±0.03 


0.142 


b 


17.8 







rl091 


16:36:49.68 


41:02:42.80 


19.404 


24.454 


0.45 ±0.10 


0.057 


d 


> 10 







rll07 


16:36:43.15 


41:02:46.79 


18.576 


23.787 


0.45 ± 0.05 


0.133 


d + p 


0.9 


0.12 




rlll4 


16:36:56.54 


41:02:47.07 


19.105 


24.265 


0.39 ±0.19 


0.150 


b + p 


0.4 


0.50 




rll27 


16:36:42.39 


41:02:54.31 


18.376 


23.831 


0.50 ±0.07 


0.270 


d 


1.3 





aso 


rl437 


16:36:53.52 


41:04:50.34 


19.283 


24.667 


0.36 ±0.06 


0.156 


d 


0.02 
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6 MORPHOLOGY AND RADII OF EROS 

6.1 Profile Fitting and Classification 

We investigate the morphology and angular sizes of the 
EROs, to as faint a limit as possible for our ground-based 
data. After some experimentation, it was found that use- 
ful morphological information could be extracted from the 
UFTI (but not Ingrid) images of galaxies to K ~ 19.5. 
Hence, this Section concentrates on the subsample of 32 
EROs on the UFTI mosaic with K < 19.5. 

Table 3 gives positions and magnitudes of these objects, 
labelled with their detection numbers in the SExtractor cat- 
alog (prefixed V for red). Using IRAF 'isophote. ellipse', a set 
of elliptical isophotes were fitted to the UFTI images of each 
ERO, using centroids, ellipticities and position angles from 
the SExtractor catalog as the starting parameters. The sur- 
face brightness (SB) on each fitted isophote, as a function 
of semi-major axis, provides a radial intensity profile. 

For each ERO, we estimate a morphological classifica- 
tion (i.e. bulge, disk or point-source) and half-light radius, 
r M , but this is complicated by the need to correct for atmo- 
spheric effects. The seeing point-spread function averaged 
0.67 arcsec FWHM, but with significant variation between 
the 16 frames of the UFTI mosaic. Hence to correct for this, 
a bright (but not saturated), relatively isolated, star was 
identified on the same frame as the each ERO. Using IRAF 
'mkobjects', a grid of model profiles was generated, includ- 
ing 

(i) a point-source, 

(ii) a set of bulge profiles, 

I(r) = I exp[-7.67( — )"*] 

(iii) disk (exponential) profiles, 

T 

I(r) = I exp - ( ) 

covering a wide range of r M in steps of 0.01 arcsec (for a 
disk profile rui = 1.679r exp ). This grid was convolved with 
the seeing point-spread function as represented by the star 
(no sky noise is added at this stage). Using 'isophote. ellipse', 
isophotal profiles were extracted from each convolved model. 
Each is normalized to the same total intensity as the ob- 
seerved ERO profile, and x 2 -tested against it. The model 
giving the smallest \ 2 is (after checking the fit by eye) then 
adopted as the estimate for morphological type and seeing- 
corrected rui- 

To estimate an error for r M , a second grid was gener- 
ated, containing multiple copies of the normalized, best-fit 
model profile. This was convolved with the seeing, and noise 
equivalent to the sky noise added. Isophotal profiles were ex- 
tracted from the resulting set of noisy models, and x 2 tested 
against the observed ERO profile, and the scatter in the re- 
sulting x 2 , °"(x 2 ) calculated. Returning to the first grid of 
noiseless models, the change A(rw) away from the best-fit 
model that increases the x 2 by a (x 2 ) ls than taken as the 
la error on r M . 

Where possible both a disk and a bulge model are fitted, 
so that the difference between their x 2 , in terms of cr(x 2 )i 
represents the significance by which one profile is favoured. 
Some EROs were not well fit by either model, but did give 
a lower x 2 when a central point-source component was in- 
cluded in the models. For these, point sources with flux frac- 
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Figure 11. Greyscale A"-band (UFTI) images of 6 of the brighter 
EROs, showing 5.45 X 5.45 arcsec areas. The galaxy r270 appears 
to be a normal elliptical, r407 a face-on disk interacting with a 
fainter galaxy, r608 a high-SB elliptical at the centre of a possible 
cluster, r622 another merging system with a double nucleus, r660 
an extended, low-SB disk galaxy, and r675 an apparently regular 
spiral. 



tions f p were combined with bulge and disk models, and we 
determined the (r/,j, f p ) combination which minimized x 2 - 
In some other EROs, the best-fit x 2 remained large because 
they have 'peculiar' morphologies, such as a double nucleus. 

For each of the 32 bright EROs, the best-fit rm, ellip- 
ticity (from SExtractor), best-fit profile type and the sig- 
nificance by which it is favoured, are given in Table 3, and 
any obvious irregularity is noted. Figure 11 shows greyscale 
images of 6 of these EROs, representing a range of morpho- 
logical types, and Figure 12 their radial intensity profiles 
with best-fit models. 



6.2 Morphology: Summary of Results 

(i) One of the 32 bright EROs, r561, appears to be a pure 
point-source, and is probably a red Galactic star. 

(ii) Of the remaining 31 EROs, 19 were best-fitted with 
bulge and 12 with exponential profiles. For some EROs the 
difference in x 2 between the two models is very small - either 
they have intermediate profiles, and/or the signal/noise is 
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insufficient for classification. Of those where one profile is 
favoured by > la, we classify 12 as bulges and 8 as disks. 

(iii) A relatively small fraction of the EROs appear to 
be ongoing mergers - we find 2/31. A further 6 (4 bulge, 2 
disk) have some visible asymmetry, most likely from recent 
interactions. 

(iv) The profiles of 9/31 EROs, including the 'cluster' 
ERO, r608, were better fitted with the addition of central 
point-source components to the underlying disk or bulge 
galaxies, although higher resolution data will be needed to 
confirm this. 

The EROs at K < 19.5 appear to consist of about a 3:2 
mixture of elliptical and spiral types, with about 1/4 show- 
ing evidence of ongoing or (more often) recent interactions. 
This is consistent with the findings of Moriondo et al. (2001) 
and Stiavelli and Treu (2001). 



6.3 Angular Sizes of EROs 

Figure 13 shows the best-fit half-light radii, tm, of the 31 
(non-stellar) EROs against K magnitude, compared to mod- 
els based on the radii of local galaxies. The size-luminosity 
relation of local E/SOs can be represented by the dou- 
ble power-law relation of Binggeli, Sandage and Tarenghi 
(1984); for M B < -20 (with h 50 = 1) 

log (r w /kpc) = -0.3(M B + 18.57) 

and for M B > -20 

log (r w /kpc) = -0.1(M fl + 15.70) 

An E/S0 galaxy with M B = -21 (M K = -25.04) at z = 
would, with the above size relation, have r M = 5.36 /igQ 1 
kpc. At z = 1, with the luminosity evolution of our passive 
model, it would have K = 18.22, and an angular ru of 0.48 
arcsec. At z = 2 it would appear fainter, K = 19.44 but 
little different in size, r M = 0.46 arcsec. Figure 13 shows the 
whole tm-Mb relation evolved to z = 1 and 2. 




r660 K=19.25 
disk r M =1.42 arcsec 



r675 K=18.70 
disk r.,-0.61 arcsec 



"..I 



Radius (arcsec) 

Figure 12. Radial (major axis) intensity profiles of the six EROs 
shown on Figure 11, compared with the best-fitting disk (dotted) 
or bulge (dashed) model profiles, and the point-spread function 
(solid line) as determined from bright stars on the same CCD 
frames. 



The brightest of these EROs, r955, has too high a sur- 
face brightness for an elliptical at z ~ 1. As it is also only 
just red enough in to have been included, the most likely ex- 
planation is that it is at a much lower redshift. However, the 
other 18 bulge-type EROs all have sizes consistent with pas- 
sively evolving E/S0 galaxies at 0.9 < z < 2. Most (12/19) 
are concentrated near the z ~ 1 locus with only two (r525, 
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Figure 13. Half-light radii rhl, as estimated by fitting radial pro- 
files, of the 32 UFTI-field EROs with K < 19.5 (symbols indicat- 
ing profile type), against K magnitude - compared with predicted 
sizes for E/SO (dashed) and disk (dotted) galaxies at z = 1 and 2, 
based on the size-luminosity relations of local E/SOs (Binggcli , 
Sandage and Tarcnghi 1984) and spirals (Cross and Driver 2002) 
with passive luminosity evolution. 

r612) on the z ~ 2 relation, as would be expected from the 
shape of N(z) (Figure 10). 

Disk EROs are more likely to be dsfEROs, and hence 
their surface brightness may be increased or decreased by 
starbursting and/or dust. Cross and Driver (2002) determine 
a bivarate brightness function for 45000 disk galaxies in the 
2dFGRS, which gives a mean effective surface brightness 
= 22.45 B mag arcsec^ 2 at L* (M B = -21.23 for h 50 = 
1), with a positive correlation between surface brightness 
and luminosity. Their best-fit relation corresponds to 

log (r w /kpc) = -0.144M S - 2.034 

with scatter it (log r M ) — 0.103. giving, for example r M = 
9.77h^ kpc for M B = -21, hence 0.87(0.83) arcsec at 
z = 1(2). A galaxy (of any morphology) in which star- 
formation is truncated long before the epoch of observation 
will be undergoing approximately the passive (E/SO) model 
luminosity evolution, so with Mb = —21 locally would have 
K = 18.22(19.44) at z = 1(2). 

Figure 13 shows the r M — Mb loci of this 'passive disk' 
model at z = 1 and 2. Most (10/12) of the disk EROs lie 
below the z = 1 locus, so if they are at z > 1 their intrinsic 
SB is greater than this model (some very dusty galaxies 
may be included as EROs with redshifts as low as z ~ 0.5 - 
see Figure 1). The two exceptions are the interacting r407, 
which lies on the 2 = 2 locus, and the very large r660, which 
is of a lower SB than this model, suggesting it is strongly 
dust-reddened (see Section 7). 

Using the PLE and M-DE models, we predict r w dis- 
tributions for bulge and disk EROs. Figure 14 compares the 



Figure 14. Angular size distribution of K < 19.5 bulge-type 
EROs in our sample (solid histogram), and Moriondo et al. (2000) 
(dotted histogram), compared to PLE (dotted), M-DE (dot-dash) 
and NE (long-dashed) models. 



N(rhi) of the K < 19.5 bulge-type EROs in our sample 
and Moriondo et al. (2000) with the models. PLE predicts 
(z) — 1.72 and (r w ) = 0.55 arcsec for the bulge EROs, wh- 
weras with M-DE they are a factor 4.07 less numerous and 
have (z) = 1.42 and (rhl) = 0.43 arcsec. With the NE model, 
the EROs are just 4 per cent less fewwr than with M-DE, 
and havewith (2) = 1.11 and (rui) = 0.50 arcsec For the 19 
bulge EROs in our sample, (ru) = 0.33 ± 0.03 arcsec, but 
the 16 of Moriondo et al. (2000) have (r M ) = 0.49±0.07 arc- 
sec (with most of the difference due to two very large cluster 
ellipticals). The combined sample would be reasonably con- 
sistent with the M-DE model, but apparently inconsistent 
with NE which overpredicts the sizes. 

Figure 15 shows model N(rhi) for disk-type EROs, 
which adopt the spiral galaxy if -band luminosity function. 
Again the evolving models assume passive luminosity evo- 
lution. As our models do not predict the fraction of disk 
EROs we have multiplied by arbitrary normalizations of 
0.15 for PLE and 0.55 for M-DE and NE. These models 
are compared to histograms of the 12 disk- type EROs in 
our sample, and the three K < 19.5 EROs from Moriondo 
et al. (2000) that were classed as disks. The models pre- 
dict (r hl ) = 0.94 arcsec for PLE, 0.76 arcsec for M-DE, 
and 0.90 arcsec for NE. In our sample the disk EROs have 
{rhl) = 0.58 ± 0.09. Again the merging model is closest, 
although these disk EROs tend to be even more compact. 

If this size discrepancy is significant, it could be the re- 
sult of size evolution - there is some evidence that z > 1 
disk galaxies in general tend to be smaller in rni than lo- 
cal counterparts of similar mass (Roche et al. 1997, 1998; 
Giallongo et al. 2000). Roche et al. (1998) hypothesised 
that this resulted from 'inside-outwards' disk formation, 
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Figure 15. Angular size distribution of K < 19.5 disk-type EROs 
in our sample (solid histogram), and Moriondo ct al. (2000) (dot- 
ted histogram); compared to PLE (dotted) and M-DE (dashed) 
models. 

with the star-formation becoming more centrally concen- 
trated at higher redshifts. Disk r M evolution appears to be 
quite moderate (< 0.1 dex) at z ~ 1, but could be suffi- 
cient for consistency with the M-DE model. It is also likely 
that many EROs classed as disks are actually disturbed, 
post-interaction galaxies which may be the process of trans- 
forming into bulge galaxies, and hence have become more 
centrally concentrated than normal spirals. 

6.4 The EROs in a possible Cluster 

In Section 3.3 we noted a possible grouping of EROs centered 
on the X-ray and radio luminous r608. There are 11 K < 21 
EROs within a 45 arcsec radius, of which 7 have K < 19.5 - 
r552, r581, r608, r626, r629, r642, and r675, an overdensity 
compared to the 1.6 K < 19.5 EROs expected by chance in 
this area. If the EROs in this area belong to a single cluster, 
they should trace an iso-redshift locus on the r\a — K plot. 
On Figure 17, the five bulge-type EROs lie neatly on the 
z = 1 E/SO locus, and the two disks are 0.1-0.2 dex below 
the z — 1 disk model. If these disks have a similar size 
offset from this model as the distribution on Figure 16, these 
radii are entirely consistent with all 7 EROs belonging to a 
single, z ~ 1 cluster. Of course, to confirm this will require 
spectroscopic redshifts (which we plan to obtain). 



7 RADIO, X-RAY AND SUB-MM 
OBSERVATIONS 

The ELAIS N2 field has also been observed: 

(i) At radio (1.4 GHz) frequencies using the VLA, 



Figure 16. Half-light radii r^i, as estimated by fitting radial 
profiles, of the 7 EROs with K < 19.5 within a 45 arcsec radius 
of r608, against K magnitude. Galaxies are labelled with detec- 
tion numbers, with symbols indicating profile type. The plotted 
models arc as in Figure 13. 



Table 4. Radio and/or X-ray fluxes of the seven K < 19.5 UFTI 
subsamplc EROs detected in the VLA and/or Chandra surveys. 



ERO 


F(1.4 GHz) (//Jy) 


F(0.5-8.0 keV) 




peak 


integral 


10~ 15 ergs s _1 cm" 2 


r256 


32.4 ±9.5 


25.1 ± 13.9 




r518 


27.7 ± 9.5 


30.0 ± 17.3 




r552 


39.3 ±9.5 


21.8 ± 11.3 




r608 


4477. ± 9.5 


5074. ± 17.8 


3.060 ± 0.891 


r642 


58.3 ±9.5 


46.7 ± 14.2 




r660 


29.1 ± 9.3 


39.3 ± 19.8 




r952 


48.0 ±9.5 


43.4 ± 15.4 





reaching a 3cr threshold for source detection F(1.4 GHz) ~ 
27.6/iJy, with resolution 1.4 arcsec (see Ivison et al. 2002). 

(ii) In X-rays using the Chandra satellite, for 75 ks 
in August 2000, reaching source detection limits F(0.5- 
8.0 keV) ~ 1.5 x 10 -15 ergs s _1 cm~ 2 with sub-arcsec reso- 
lution ( Manners et al. 2002). 

(iii) At 850pim using SCUBA, reaching a 3.5ct source 
detection limit of F(850/im) ~ 8 mjy (Scott et al. 2002). 

Here we discuss these observations for the K < 19.5 
UFTI EROs only. Seven of these 31 galaxies are detected at 
1.4 GHz sources, only one as an X-ray source, and none in 
the sub-mm. Table 4 gives the fluxes of the detections. In 
r256, r518, r552, r642 and r952, the radio emission is rela- 
tively weak and consistent with a point source, and the host 
galaxies appear to be regular ellipticals (the outer regions of 
r518 are mildly disturbed). The radio emission is stronger 
than would be expected for typical L* ellipticals (< lO^iJy 
at all z > 1) and seems more likely to be produced by weak 
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Figure 17. Radio contour map centred on r660, generated from 
VLA data using the AIPS package, superimposed on the UFTI 
A"-band image. The radio coutours are at 2 and 3a where a = 
9.2/iJy/beam (dotted contours are negative). 



AGN, below the Chandra threshold, rather than powerful 
starbursts. 

The galaxy r660, similar to these five in radio flux, dif- 
fers in that (i) the radio emission is elongated, with major 
axis FWHM ~ 2.3 arscec on a position angle 27±18 degrees, 
and ( Figure 18) aligned with the K-b&nd long axis, (ii) the 
galaxy is a large low-surface brightness disk. This suggests 
that the radio emission is produced by an extended dust- 
reddened starburst. With the Condon (1992) radio SED, the 
39.3/iJy flux corresponds to a rest- frame luminosity (yL v ) at 
z = 1(2) of L1.4 = 10 39 - 66 (10 40 - 34 )^5 2 ergs s" 1 . If L1.4/SFR 
is in the range bracketed by the Carilli (2000) and Con- 
don (1992) relations, as in Roche et al. (2002), the total 
SFR = 99-233(474-1115)ft5 2 Moyr" 1 at z = 1(2), and r660 
is a very powerful starburst galaxy, probably similar to local 
ULIRGs. 

The non-detection of any of these 31 EROs in the 
SCUBA survey sets an upper limit on their SFR. Scott et al. 
(2002) estimate the detection limit of F(850/im) ~ 8mJy to 
correspond to a SFR ~ 1800/ig, 2 Moyr" 1 , at all 1 < 2 < 10, 
for the Salpeter IMF, which for the IMF of our models be- 
comes ~ 1200ft5 2 M yr _1 . Most of these bright EROs, in- 
cluding the double nucleus and disturbed galaxies, are not 
detected in the VLA data either, which sets stricter upper 
limits on SFRs, ~ 70-164(335-785)/i^ 2 Moyr" 1 at z = 1(2). 
As most of the brightest EROs will be closer to z ~ 1, this 
implies that the average SFR of dsfEROs at these redshifts 
is less than ~ 200/i^ 2 MQyr -1 . Nevertheless, a significant 
minority - e.g. r660, and the numerous examples of fainter 
EROs associated with SCUBA sources (Smail et al. 1999; 
Dey et al. 1999; Ivison et al. 2001; Lutz et al. 2001) - are 
real ULIRGs. 

By far the strongest radio source in this sample, at 
5m Jy, is r608, the X-ray source and possible cluster central 
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Figure 18. Radio contour map centred on r608, superimposed 
on the UFTI X-band image. The radio contours are at 2,3,4...10<x 
where a = 9.2/iJy/bcam. 

galaxy. The X-ray and VLA data (Figure 19) show most of 
the emission from a point-source concentric with the if-band 
profile. The radio flux is far too strong to be attributed to 
star-formation (the non-detection with SCUBA limits the 
starburst contribution to < 1 per cent of this). If r623 is 
at z ~ 1, the rest-frame radio luminosity, L1.4 — 10 41 ' 75 
ergs s" 1 . This is 3 orders of magnitude above that of normal 
spirals but 3 orders below the most powerful radio galax- 
ies. If we assume passive optical evolution and (1 + z) 3 ra- 
dio evolution for this galaxy then at z = it would have 
M R = -23.24 and L 1A ~ 10 40 ' 85 ergs s" 1 , or P 1A = 10 24 ' 70 
W Hz -1 , and its position on the radio-optical luminosity 
plane (Ledlow and Owen 1996) would then be typical of 
local FRI radio galaxies. 

The X-ray flux, F(0.5-8.0 keV) = 3.06 x 10~ 15 ergs 
s _1 cm -2 , corresponds at z ~ 1, assuming a /„ oc spec- 
trum, to L x (0.5-8.0keV) = 10 43M ergs s" 1 , which could 
also be expressed in the vL v form as LikcV = 10 43 ' 27 
ergs s~ x . From this, the radio to X-ray spectral index, 
a RX ~ 41.75-43.27+8.24 = °- 815 > which is also consistent 
with an FRI radio galaxy (e.g. Capetti et al. 2002). 



8 SUMMARY AND DISCUSSION 
8.1 Summary 

(i) The number counts of EROs are significantly lower than 
expected if all local E/SOs evolved as in a PLE model, but 
are much higher than the predictions of CDM-based hierar- 
chical merging models (Cimatti 2002b). We find (in agree- 
ment with Firth et al. 2002) that the ERO counts are consis- 
tent with completely non-evolving model for E/SOs. We can 
also fit the ERO counts with a physically plausible evolving 
model, 'M-DE', in which both the characteristic mass and 
the comoving number density of passively evolving galax- 
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ies decrease with redshift. Physically this evolution could be 
accounted for by the continual formation of new passive el- 
lipticals from the mergers and interactions of star-forming 
spirals. 

Our best-fitting 'M-DE' model incorporates an 
observationally-estimated merger rate (Patton et al. 2001), 
with the negative density evolution parameterized by R^ — 
—0.46, meaning that only ~ 35 per cent of the present-day 
comoving number density of E/SO galaxies formed at z > 3. 
This fraction is higher than the ~ 15 per cent estimated by 
Stiavalli and Treu (2001), because we have (a) included all 
types of ERO as possible E/S0 progenitors (b) included a 
reduction in mean mass galaxy and hence L* at high red- 
shift, due to merging thus allowing a higher number density 
($*)of passive galaxies. 

(ii) We examined the morphology of a bright (K < 19.5) 
subsample of 32 EROs on the UFTI mosaic. One appeared 
to be a Galactic star, the other 31 a ~3:2 mixture of elliptical 
and disk profile types. About 1/4 showed some evidence 
of asymmetry or disturbance, e.g from recent interactions, 
but only two were obvious mergers with double nuclei. This 
mixture of morphologies is consistent with the findings of 
Moriondo et al. (2001) and Stiavelli and Treu (2001) for 
similar samples of EROs. 

By fitting radial profiles, we estimated seeing-corrected 
half-light radii. The r M of the bulge-type EROs were con- 
sistent with passively evolving E/S0 galaxies in the redshift 
range (0.9 < z < 2.0) predicted by the M-DE model. This 
is in agreement with Moriondo et al. (2000) who found the 
surface brightness-size relation of six bulge-type EROs at 
known rcdshifts of z ~ 1.3 to be consistent with passively 
evolved (A(B) — —1.4 mag) ellipticals. 

Although both the M-DE and non-evolving models fit 
the ERO counts, these models give different angular size 
distributions, and our observed ERO radii are much more 
consistent with the smaller radii predicted by M-DE. This 
implies that, whatever the L* evolution, the surface bright- 
ness evolution of EROs is at least that predicted for passive 
evolution, and if the L* evolution is much weaker this must 
be the result of merging. The average r^i of the disk-profile 
EROs is even smaller than the M-DE prediction and possible 
explanations include inside-outwards disk formation and/or 
morphological evolution. 

(iii) Radio observations with the VLA detect emission above 
the 3<r limit of F(1.4 GHz) ~ 27.6^Jy for 7 of the 31 ERO 
galaxies with K < 19.5 on the UFTI field. The strongest 
of these sources, at 5m Jy, would be a radio galaxy, perhaps 
like local FRIs if these evolve as strongly as L ra d oc (1 + z) A . 
The host galaxy is also a Chandra X-ray source and may 
be the central giant elliptical of a cluster of EROs. Of the 
other, much fainter (30-60/iJy) detections, five are point 
sources within apparently regular galaxies, and may be weak 
AGN. The seventh, the galaxy r660, shows elongated radio 
emission aligned with the K-h&nA image. The galaxy is a 
large disk of lower surface brightness than the other bright 
EROs. These properties suggest it is a powerful but very 
dust-reddened starburst, with a SFR 100-1000 M yr _1 . 

The non-detection of the other EROs, which include 
double- nucleus and disturbed objects, implies their SFRs 
are lower. This was interpreted as indicating the mean SFR 



of the dsfEROs to be < 200h^ M yr _1 . This is probably 
consistent with the Cimatti et al. (2002a) estimate of ~ 
lOO/i^ 2 M Q yr _1 , based on the mean [OIIJ3727 and UV fluxes 
of spectroscopically observed dsfEROs (see 8.2 below), after 
correcting for E(B — V) = 0.5 mag extinction. 

(iv) We investigated the angular correlation function ui(6) of 
EROs on the UFTI and Ingrid fields. Positive clustering is 
detected at > 2a for the EROs to K = 19.5-20.0, with an es- 
timated u(6) amplitude about an order of magnitude higher 
than that of all galaxies to the same K mag limits. Our ERO 
clustering results are consistent with those of Daddi et al. 
(2000) and Firth et al. (2002), and using a combination of 
these, we estimate that the brighter (K < 20.0) EROs have 
a comoving correlation radius ro — 10-13 h^Q Mpc, depend- 
ing on the strength of V evolution. Statistical uncertainties 
remain large and it there is clearly a need for further ERO 
clustering analyses with both larger and deeper samples. 



8.2 Discussion 

8.2.1 Star-forming EROs 

Cimatti et al. (2002a) , using VLT spectroscopy, classified the 
spectra of about two-thirds of a sample of bright (K < 19.2) 
EROs, and found approximately half (50 ± 17 per cent) to 
be dusty star- forming and half to be old (> 3 Gyr) pas- 
sive galaxies. The averaged spectrum of the dsfEROs showed 
[Oil] 3727 emission and Balmer absorption lines and closely 
resembled an e(a) type (Poggianti and Wu 2000) 'very lu- 
minous infra-red galaxy' with stellar reddening E(B — V) ~ 
0.5. 

The VLA data indicate that only a small fraction of 
EROs can be ULIRGs, and furthermore the proportion of 
EROs which are ongoing mergers or very disturbed is con- 
siderably lower than 50 per cent, so some dsfEROs must 
have more regular bulge or intermediate profiles. This sug- 
gests that many of the dsfEROs are in late post-merger 
stages, with low and declining SFR and increasingly reg- 
ular morphology. These galaxies may subsequently become 
E/SOs and contribute to the increase with time in the co- 
moving number density of passive galaxies. The dsfEROs 
and pEROs would then be separate stages of an evolution- 
ary sequence. At z ~ 1-2 we see the second (pERO) stage 
for the earliest formed ellipticals, contemporaneous with the 
dsfERO stage of those forming later, while deeper ERO sur- 
veys would reach the dusty starburst phase of the first mas- 
sive ellipticals, at z > 3 (i.e. the SCUBA sources). 

On the basis of the Cimatti et al. (2002a) spectroscopy, 
dsfEROs would make up ~ (0.5 ± 0.17) x 14 = 7 ± 2 per 
cent of all K < 20 galaxies and to this limit have a surface 
density 0.69 ±0.23 arcmin -2 . These high numbers imply the 
dsfERO phase is prolonged. For an order-of-magnitude esti- 
mate of the total dsfERO lifetime, tdsf, we assume that (i) 
each present-day E/S0 galaxy is associated with an average 
of one dsfERO 'event' at higher redshift, (ii) that the K- 
band magnitudes of dsfEROs can be approximated by the 
passive model with merging (as they seem to be similar in 
surface brightness to the pEROs). 

If tdsf were to last for the entire 'ERO epoch', z > 
0.93, the dsfERO count could then be represented by the 
z > 0.93 E/SOs in a model with merging mass evolution 
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(R m = 0.3) but a comoving number density remaining at 
the presentOday value (i?^ = 0). This would be intermediate 
between the merging and M-DE models on Figure 6 and to 
K — 20 gives 2.9 arcmin -2 . The observed dsfERO count is 
lower, implying a shorter tdsf. As most K < 20 EROs will 
be at 0.93 < z < 2, or lookback times in the period (10.4- 
UA)h^ Gyr, hence t dsf ~ %f x (14.4-10.4) = 1.2 h^ 1 
Gyr. 

8.2.2 Clustering and evolution of the EROs 

We estimate the intrinsic clustering of bright (K < 20) 
EROs as ro — 10-13 hio Mpc comoving, depending on 
the strength of L* evolution. This is even stronger than the 
clustering of present-day giant ellipticals, ro — 8 /i^oo, and 
implies that if the z ~ 1-2 EROs evolve into the z ~ 
E/SOs, the increase in their comoving number density must 
involve the assimilation of less clustered types of galaxy into 
the class of passive galaxies, reducing the overall ro. 

It may be significant that Abell galaxy clusters, with 
ro — 21/ijqq Mpc (Abadi, Lambas and Muriel 1998), and 
possibly SCUBA sources (Almaini et al. 2002), may be even 
more clustered than the z ~ 1-2 EROs. On the other hand, 
the Lyman break galaxies at z ~ 3 are on the whole more 
moderately clustered, with a comoving ro — 2-5 h^ Q Mpc 
(Giavalisco and Dickinson 2001; Arnouts et al. 2002), which 
suggests they are mostly the progenitors of disk galaxies. 

Daddi et al. (2002) found, on the basis of a small sam- 
ple of redshifts, some evidence that the strong ERO cluster- 
ing is associated with the pERO types only. They estimate 
ro < 2.5 /Ij'qq Mpc for dsfEROs, and argue that this may be 
evidence against dsfEROs being E/S0 progenitors. 

However, a large difference in pERO and dsfERO clus- 
tering would be expected if (i) the pEROs form at very 
dense mass concentrations in the early Universe, the 'seeds' 
of Abell clusters and SCUBA sources (and so have simi- 
lar clustering), and (ii) the addition with time of dsfEROs 
(formed in mergers of the relatively weakly clustered disk 
galaxies) to the ERO class produces a continual dilution of 
the clustering - until, by z ~ 0, the passive galaxy ro is as 
observed for local E/SOs. 

Secondly, IRAS (60^im) selected galaxies, which may be 
considered as low redshift counterparts of the dsfEROs, are 
weakly clustered (ro — 4 Ii~[q Mpc), but nevertheless have 
many properties suggesting they are destined to evolve into 
E/SOs (Genzel et al. 2001). 

Thirdly, many low-redshift E/S0 are observed to have 
have 'boxy' or 'disky' isophotes, and this may be 'fossil' 
evidence that they evolved from both types of ERO. The 
'boxy' ellipticals tend to be more massive and more con- 
fined to the centres of clusters, with slightly older stellar 
populations. The 'disky' galaxies tend to have higher rota- 
tion velocities (v ro t) but lower internal velocity dispersions 
(<Ti). These differences may result from the 'disky' ellipticals 
having formed at a relatively late epoch, when large spirals 
were already present - perhaps (at one extreme) from the 
merger of a single pair of spirals, which is likely impart a 
substantial angular momentum to the post-merger. At the 
other extreme, the oldest ellipticals might be formed through 
multiple mergers of very many small disks in the core of a 
proto-cluster, giving a summed angular momentum vector 
much closer to zero, and 'boxy' isophotes This latter process 



may have been observed in deep K images of high-z radio 
galaxies, which show massive ellipticals at z < 3 but a large 
number of clustered sub-components at z > 3 (van Breugel 
et al. 1999). 

At z ~ 1-2, the progenitors of today's 'boxy' and 'disky' 
E/SOs would then be, respectively, the pEROs and dsfEROs. 

8.2.3 Future plans 

We hope to investigate these scenarios further through 

(i) Analysing the clustering and morphologies of of 
larger samples of EROs, with multi-colour imaging (e.g. R, 
I, J, H and K) to enable separation of the dsfEROs and 
pEROs, to fainter (K > 22) limits. We now have some data 
for this. It may prove especially useful to compare the clus- 
tering properties of EROs and SCUBA sources. 

(ii) High-resolution spectroscopy of a diverse sample of 
EROs, with the aim of not only determining the N(z), but 
also of estimating properties such as current SFR, internal 
dust extinction, metallicity and kinematics, and applying 
new techniques of age-dating the stellar populations. From 
the velocities er; and v ro t, measurable from ground-based 
spectroscopy, on a v ro t-o~i plot (Genzel et al. 2001), the ERO 
progenitors of 'disky' and 'boxy' ellipticals might be identi- 
fied. 
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